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Abstract 

We propose a simple Ansatz for the three generation neutrino mass matrix 
Mu which is motivated from an SO(IO) grand unified theory. The Ansatz 
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themselves and to reconstruct, with some assumptions, the matrix Mi, . 
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I. INTRODUCTION 



In the last few years there has been another wave of excitement regarding the question of 
neutrino masses. This is largely due to the many new experiments testing neutrino oscilla- 
tions, most notably the positive indications obtained by Super Kamiokande on atmospheric 
neutrino oscillations [Q. Similar indications come from other experiments 0-^. The solar 
neutrino experiments have for many years provided independent evidence for neutrino os- 
cillations ||6|-[T0|. Accelerator and reactor experiments have also played an important role. 



They have furnished strict bounds on neutrino oscillation parameters [|TT]-[T^. In the case of 
the LSND experiment ||TB[ at Los Alamos evidence for — > Vf. oscillation has been reported. 
See refs [17| for recent reviews. 



It is hoped that new experimental results can be used to determine the neutrino squared 
mass differences and mixing angles. In turn, these may help to infer the neutrino mass 
matrix. This is presumably a possible gateway to a more fundamental theory beyond the 
standard model. Of course this is a highly speculative area, and even though there are 
many imaginative proposals |jl8|, it seems fair to say that the the true answer is essentially 



unknown. In order to make progress in this direction, it seems useful to investigate various 
plausible Ansatze for the neutrino mass matrix. From this point of view we propose the 
Ansatz for the 3 generation neutrino mass matrix, M^: 

Tr(M,) = (1.1) 

and investigate its consequences. We are considering the neutrinos to be represented by 
2-component spinors so that, in the most general situation. My is an arbitrary symmetric 
complex matrix. 

As we will see in section II, Eq. (1.1) can be motivated from an SO(IO) grand unified 
model []T9|, in which it may be derived with some assumptions. Physically, Eq. (1.1) 



corresponds to the well known approximate signature of grand unification that ~ 3. 
Furthermore we will see in sections IV and V that Eq. (1.1) can be straightforwardly 
combined with experimental information to get an idea of the neutrino masses themselves 
as well as the "texture" of M^. Relevant matters of notation are discussed in section III 
while a summary is presented in section VI. 



II. PLAUSIBILITY OF THE ANSATZ 

In the SO (10) grand unification model each generation contains one light massive two 
component neutrino and also a very heavy one which is "integrated out" according to the 
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"seesaw mechanism" p^. The effective 3x3 neutrino mass matrix takes the form: 

M^ = Ml- MIMhMd (2.1) 

where M^,, and M/j are respectively the mass matrices of the hght neutrinos, heavy 
neutrinos and heavy-hght mixing (or "Dirac matrix"). Generally the second, seesaw, term 
is considered to dominate. Here however we shall assume the first term to be the dominant 
one. This is necessary for the present derivation of Eq. (1.1) to hold. Also, a rough order of 
magnitude estimate for the second term would be ^^^SQev about 3 x 10"'^ eV. Thus, the 
seesaw term could be negligible if neutrino masses turn out to be appreciably larger than 
this value. Now in SO(IO), Higgs mesons belonging to the 10, 120 and 126 representations 



can contribute to the fermion masses at tree level. One has ||21| for the down quark, charged 
lepton and light neutrino mass matrices, 

Md = aS{10) + 6A(120) - ^cS{126) 
rM^ = aS{10) + dA{120) + cS(126) (2.2) 
sMl = eS{126) 

where a, b, c, d, e are numbers representing Higgs meson vacuum values. S(10), A(120) and 
S(126) are the matrices of the Yukawa type constants which couple the fermions to the 10, 
120 and 126 Higgs mesons respectively; the matrices S(10) and S(126) must be symmetric 
while A(120) is antisymmetric. Finally, r ~ 3 is a renormalization factor for comparing the 
quark masses with the charged lepton masses at a low energy scale rather than at the grand 
unified scale; s is a similar factor for the neutrino masses. With the stated assumption that 
the Ml term dominates in Eq. (2.1) we get 

Tr(M^) oc Tr(Md) - rTr(Me) (2.3) 

which clearly also holds when any number of lO's or 120's are present but only a single 126. 
The matrices appearing in Eq. (2.3) are so far essentially unrestricted complex ones. To 
proceed, we make the further assumption that the matrices are hermitian. Then M^ and 
Mu may each be brought to diagonal form by unitary transformations. Thus the right hand 
side of Eq. (2.3) may be evaluated to yield approximately, 

Tr(M^) oc m{b) - rm(r) ~ (2.4) 

according to a well known numerical success, based on the observation that r ^ 3, of grand 



unification |22]. Note that we have not needed to assume that the mass matrix has any 



2 



zero elements.^ Even if the cancellation on the right hand side of Eq. (2.4) is not perfect, 
it should still be a good approximation. In an SO(IO) model where the mass matrices are 
hermitian, will be real symmetric. We will investigate this case and also the possibility 
that the more general case holds. 

III. SOME NOTATION 

Our plan is to combine the Ansatz Eq. (1.1) with experimentally obtained results on 
neutrino oscillations in order to learn more about My itself. For this purpose it may be 



helpful to set down our notation for the pieces of the effective SU(2)i x U(l) theory 



involving neutrinos and to make some related remarks. 

The free Lagrangian containing three two component massive fields is: 

Cfree = -ip^O-^d^p - ^(p'^(T2M^p + h.c), (3.1) 

where = Mj is the (not yet diagonalized) neutrino mass matrix of the underlying theory 
to be identified with the matrix in Eq. (1.1). Note that we are free to multiply the first 
mass term in Eq. ( |3.1| ) by an overall arbitrary phase which is a matter of convention. It is 



possible PB| to find a unitary matrix U which brings to real, positive, diagonal form in 



the following way: 

U^MM = M, M = diag(mi, ma, mg). (3.2) 
The mass diagonal fields i> are then 

p = Uu. (3.3) 

Similarly, the column vector of left handed negatively charged leptons in the underlying 
theory, is related to the mass diagonal fields e^, by 

El = neL, (3.4) 

where fi^ = Q~^. 



*In 1 24 1 a similar mechanism was studied for Mh where, in addition, a special combined Fritzsch- 
Stech Ansatz was used. Here we are not making any special Ansatz of this type for the mass 
matrices. 
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Combining factors from Eq. (3.3) and Eq. (3.4) we obtain the unitary mixing matrix, 
K for the charged current weak interaction, 



K = n^u. 



(3.5) 



This appears in the Lagrangian term, 



Cint = -y^W'eLj^Ku + h.c. 



(3.6) 



where a conventional four component Dirac notation with 75 diagonal is being employed and 
u has only the first two components non zero. Next we parameterize K p3[. It is possible 
to restrict deti^' = 1 by adjusting an overall phase which can be absorbed in ei- Then we 
write 



where 



K = LL!Q{a)u23{l]23, 023)t^l2('7l2, 012)^13(^?13, ^is), 



cuo(«) = diag(e^"^e^"^e^°^), 



(3.7) 



(3.8) 



with 0:3 = — («! + 0:2) and, for example 



'^12(^712, 012) 



cosr7i2 e*''^^^sinr7i2 



p-«<Pl2 



smr7i2 







cosr7i2 
1 



(3.9) 



Eq. (3.7) contains the eight parameters needed to characterize an arbitrary unitary uni- 
modular matrix. From the standpoint of Eq. (3.6) it can be further simplified by using the 
freedom to rephase eLUjQ^{a) without changing the free part of the charged lepton 

Lagrangian. On the other hand, the form of the mass terms in Eq. (|3.1| ) shows that the 
neutrino fields can not be rephased. Thus a suitable minimal parameterization]] for K in 
(3.6) is 



K = UJ23{V23, </'23)^12('7l2, 012)^^13 (^713 , 013), 

involving three "angles", rjab and three "phases", (j)ab- Note the identity 

UJo{a)uJab{'nab, (l)ab)i^o^{a) = ^abijlab, Ola + 0afe " Oib) ■ 



(3.10) 



(3.11) 



^This is written out in detail in Eq (2) of 
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This identity may be used to transfer two of the phases (pab in Eq. (3.7) to a diagonal matrix 
on the right of K as, for example, 



K = uj23{r]23, 0)u;i2(?7i2, 0)u;i3(?7i3, 0i3)c^o ^(^)' (3-12) 

where ti + T2 + = 0, which may be used instead of Eq. (3.10). 

We also need the formula for the amplitude of neutrino oscillation. For the case when a 
neutrino, produced by a charged lepton of type a, "oscillates" to make at time t, a charged 
lepton of type b, we have 

amp(a ^ 6) = E KaKbaC-'"'"', (3.13) 

a 

where the sum goes over the neutrinos of definite mass, rria- Inserting the parameterization 
Eq. (3.12) into Eq. (3.13) shows that the effect of the factor uJq'^{t) cancels out. Thus for 
ordinary oscillations, K is parameterized by three angles and one CP violating phase as for 
the CKM quark mixing matrix. On the other hand, the two additional CP violating phases 
Ti and T2 show up if one considers neutrino-antineutrino oscillations or neutrinoless 
double beta decay The formula for the probability, Pab is gotten by taking the squared 
magnitude of Eq. (3.13) and replacing the exponential factor Eat by {E + 171^/ {2E))L, 
where E is the neutrino energy and L is the oscillation distance. For practical reasons it is 
very important to take account of the experimental uncertainties in E and L. The simplest 
approximation [PS| is to define b = L/{4E) and assume that one can smear Pab with a 



Gaussian distribution in b. bo is defined as the mean value and cr^ as the standard deviation 
appropriate to the particular physical setup. Then we find for the smeared probability 

{Pab) = 5a6 - 2 E [Re{fapab)il - cos(26om5„)exp(-2a,' (mj,)') ) (3.14) 

Q</3 

+ Im ( fafSab ) sm{2boml ) exp ( - 2o-^ {ml ) ^ ) ] , 



where fa^ab = KaaKlf^Kl^Kbp and = m| - m^. Notice that when Im(/„^ab) = 0, {Pab) 



is independent of the sign of m\ 



IV. LEARNING ABOUT FROM EXPERIMENT 

Since Tr {My) = provides only two real equations for 12 real parameters, it is clear that 
it has a relatively small amount of predictivity. In particular it cannot say much about the 
texture (e.g. possible zeroes) of which is suppposed to derive from a deeper theory than 
the standard model. On the other hand, we shall see that our Ansatz is complementary to 
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the results which should emerge from analysis of neutrino oscillation experiments. Together, 
they should enable us to actually (with some conditions) reconstruct M^^. 

First, in this section, we shall consider My to be hermitian so that the argument in favor 
of Tr {My) = presented in section II holds without any further assumptions. Since My is 
symmetric it must be real. It can be brought to real diagonal form via a real rotation R 
as R^MyR. However there is no guarantee that all eigenvalues of My will be positive. We 
can make them all positive by rephasing the diagonal fields with negative eigenvalues by a 
factor i (see Eqs. (3.1) and (3.2)). This means that the general diagonalizing matrix U in 
Eq. (3.2) now takes the form 

U = Ruj, (4.1) 

11 1 
where Uab = SabVb with rj^ = 1 for a positive eigenvalue and 1]^ = i for a negative eigenvalue. 

We notice that Eq. (4.1) is of the form Eq. (3.12) for which we already noticed that the 

factor u cancels out in the neutrino oscillation formula Eq. (3.14). Furthermore only the 

square of the mass is relevant in Eq. (3.14). Thus we choose to work in this section with 

some negative masses and no factor u in Eq. (4.1). 

To avoid confusion, we remark that the factor u in Eq. (4.1) does not introduce any CP 

violation in the theory ||2^ since My is real in any event. 



Now let us suppose that an experimental analysis of all neutrino oscillation experiments is 
made based on a formula like Eq. (3.14) (or one which treats the experimental uncertainties 
in a more sophisticated way). Furthermore assume that the CP violating phase 0i3 in Eq. 
(3.12) is negligible. Then we should know the magnitudes of the squared neutrino mass 
differences 

{m^Y - (mi)' = A, (ma)' - (ms)' = B, (4.2) 

where A and B can be either positive or negative. Then, assuming the leptonic theory to 
be CP conserving, our Ansatz would imply 

= Tr (My) = Tr (RMR^) = mi + ms + mg, (4.3) 



where Eq. (3.2) was used. Eqs. (4.2) and (4.3) comprise three equations for the three 
neutrino masses mi, m2 and ms. We can solve to get: 

B^ 



Va^TWTab -(f^ ■ 2 

B - (mi? 
2mi 

m3 = —mi — m2- 



(4.4) 
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This leads to a limited number of solutions, depending on sign choices. 

If we make the further assumption that the charged lepton mixing matrix Q in Eq. 

(3.4) is approximately the unit matrix (this is expected to be a reasonable but not perfect 
approximation) we can identify R ^ K which would be obtained from experiment. Then, 
using the masses found in Eq. (4.4), we could reconstruct as 

^ RMR^. (4.5) 

To proceed, we need only insert the experimental results for A, B and K in (4.4) and 

(4.5) . Of course, it is presumably the task of the next decade to solidify the experimental 
determination of these quantities. We can, at the moment, only give a preliminary discussion. 
For this purpose we will use the results of a recent preliminary analysis of all neutrino 



experiments by Ohlsson and Snellman ||3^. These authors found, by a least square analysis, 
a best fit for (our notation) 1^41, \B\ and the leptonic mixing matrix K. They used the 
formula Eq. (3.14) with a suitable choice of 6o and cr^ for each experiment. Furthermore, 
they made the simplifying assumption that K is real. Finally they only searched for a fit in 
the range 10"^ eV^ < \A\ < lO'^ eV^ 0.2 eV^ < |S| < 2 eVl This range corresponds to 
mass difference choices for which the MSW effect |^ for solar and atmospheric neutrinos 
is not expected to be important and so greatly simplifies the analysis. Thus there is no 



guarantee that the solution of |^0| is unique. Altogether they fit sixteen different solar 



neutrino, atmospheric neutrino, accelerator and reactor experiments, including LSND. The 
best fit is: 

\A\ = 2.87 X 10"^ eV^ \B\ = 1.11 eV^ (4.6) 
for the squared mass differences and 



K 

^ ^ exp 



0.7052 0.7052 0.0732 
-0.6441 0.5940 0.4820 
0.2964 -0.3871 0.8731 



(4.7) 



for the lepton mixing matrix K. As discussed above we will identify K ^ R here, keeping 
K real but allowing for negative masses. The best fit matrix K was obtained to be similar 
but not identical to the "bimaximal mixing" matrix ||32|| . 

With the best fit squared mass differences in Eq. (4.6), our model predicts, from the 
first of Eq. (4.4), eight different possibilities. These correspond to four different sign config- 
urations for A and B times the two possible signs for mi. However, only two of these eight 
are essentially different; these are 
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type I : mi = 0.6082 eV, ma = 0.6084 eV, mg = 
type II : mi = 1.053701 eV, ma = -1.053565 eV, 



-1.2166 eV (4.8) 
= -0.000136 eV (4.9) 



The other solutions correspond to interchanging whichever of |mi| and |m2| is greater (which 
has only a negligible effect since they are almost degenerate) or reversing the signs of all 
masses. Physically it is clear what is happening: the smallness of |y4| compared to \B\ in 
Eq. (4.2) forces |mi| ^ |m,2|. Then we have either mi ^ m2 with, using the constraint Eq. 
(4.3), ms ^ — 2mi or mi ^ — m2 with very small. 

Since we have assumed the neutrinos to be of Majorana type for our plausibility argu- 
ment in section II, their interactions will violate lepton number. Then they should mediate 
neutrinoless double beta decay (PPou) Such a process has not yet been observed and 
an upper bound has been set for the relevant quantity 

3 



la rric 



a=l 



(4.10) 



The best upper bound at present is ||3^ (m^) < 0.2 — 0.6 eV, reflecting some uncertainty in 
the estimation of the needed nuclear matrix elements. 

Substituting the best fit for the matrix K from Eq. (4.7) together with our results in 
Eqs. (4.8) and (4.9) into Eq. (4.10) yields predictions for the two cases: 



type I : (m^) = 0.60 eV, 
type II : (m^) = 6.7 x 10"^ eV. 



(4.11) 
(4.12) 



Both solutions seem to be acceptable, the type I case marginally but the type II case 
definitely. Note that the small value for (m,^) in the type II case is due to the best fit 
prediction Ku = K12 and also to the fact that m2 is negative. The same value would 
clearly result if we made m2 positive and set i^i2 = 0.7052i as discussed around Eq. (4.1) 
above. 

Finally, let us reconstruct the underlying neutrino mass matrices for each of the two 
cases. We use Eq. (4.5) based on the assumption that My is real and also our ansatz to find 
(in units of eV): 



type I : = Ml 



type II : 



M, 



T 



0.5985 -0.0643 -0.1167 
-0.0643 0.1843 -0.7680 
-0.1167 -0.7680 -0.7828 

6.7 X 10-5 -0.9199 0.5078 
-0.9199 0.0654 0.0410 
0.5078 0.0410 -0.0654 



(4.13) 



(4.14) 
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The type I matrix does not have an excellent candidate for a "texture" zero. However the 
small value of {M^)-^^ in the type II case is certainly suggestive. These matrices lead to 
neutrino masses and a mixing matrix which give a best fit to all present data. It will be 
interesting to see if either of them hold up in the future. 

Incidentally, on comparing Eqs. (4.13) and (4.14) it is amusing to observe the large 
difference in two mass matrices "generated" in the same way except with respect to how 
Tr (M^) = is satisfied. 

V. CASE OF COMPLEX 

It seems interesting to also investigate the Ansatz Tr{M,^) = when is no longer 
restricted to be real. This also raises the problem of constructing the unitary diagonaliz- 
ing matrix U in Eq.(3.2), in terms of the experimentally measured lepton mixing matrix 
Kexp- For simplicity, as before, wc will make the approximation that the charged lepton 
diagonalizing matrix Q is the unit matrix. 

Apart from an overall (conventional) phase we may write 

U ^ uJo{a)KexpUJo^{r), (5.1) 

where the 2-parameter quantity cuq was defined in Eq. (3.8). Since Kexp has four parameters 
U in Eq. (5.1) is described by eight parameters. As mentioned before, the two parameters 
in u!q^{t) are not measurable in neutrino oscillation experiments but show up when one con- 
siders {PPi)ov The two parameters in ujq{(t) may be eliminated, for experimental purposes, 
by rephasing the charged leptons. However, for the theoretical purpose of reconstructing 
the underlying neutrino mass matrix M^,, their existence cannot be ruled out. (They also 
do not contribute to {PP)oi,.) 

For the purpose of relating the Ansatz on M^, to the physical neutrino masses in M, we 
note 

IV(M,) = IV(X->o"'(2<7)i^eVo(2T)M). (5.2) 

For further simplicity of the analysis we adopt the special case u!Q{2a) — 1 and also 
identify K^xp with the real best fit in Eq. (4.7); our Ansatz now reads 

= Tr(M^) = Tr [ujo{2t)m) . (5.3) 

With the redefinitions Pi — 4ti + 2t2 and P2 — 2ri + 4t2, Eq. (5.3) becomes 
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FIG. 1. Geometrical picture of Eq. (5.4). 







(5.4) 



This may be conveniently visualized as the vector triangle shown in Fig. |1]. 

Combining Eqs. (5.4) and (4.2) gives four real equations for the five unknown quantities 
(mi, m2, ms, /?2)- Thus we have (for each set of {A,B) sign choices) a one parameter 
family of solutions. It is convenient to choose this parameter to be m^. Then rrii and m2 
may be found from the equations (4.2), provided that solutions exist. In this way all three 
sides of the triangle in Fig.|l| are determined. The angles may finally be found as 



COS/?2 



TTln 



2m2m3 



smpi = smp2. 

mi 



(5.5) 



We also need to investigate the constraint arising from the non-observation of (/?/3)oy 
Eq. (4.10) now becomes, with Eq. (5.1) as the mixing matrix 



(m^) = mi{Kexp ii) e + m2{Kexp 12) e + m^i^Kexp 13) 
Using the Ansatz constraint Eq. (5.4), Eq. (5.6) may be rewritten as 



{m^) = {Kexp 12)^ - {Kexp nf 



exp 13 ) 



exp 11 y 



1713 



(5.6) 



(5.7) 



This form is very convenient when identifying K^xp with the best fit solution in Eq. (4.7). In 
the present context such an identification corresponds to CP violation for the {PP)ou process 
but not for usual neutrino oscillations. Since the (11) and (12) matrix elements are equal in 
Eq. (4.7) we find the simple result 
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(m^) = 0.49m3. (5.8) 

Thus if the upper bound on {niy) is conservatively identified as in the 0.2 — 0.6 eV range, 
we should have in this case 

mg < 0.41 - 1.22eV. (5.9) 

In the present complex case there is a continuum of possible solutions labelled by those 
values of satisfying Eq. (5.9), rather than just the two possibilities found in Eq. (4.8) 
and Eq. (4.9). Actually, the continuum separates roughly into two classes similar to either 
Eq. (4.8) or Eq. (4.9). In the generalized type I class, is of the order B = 1.11 eV while 
mi and m2 are related to nearly oppositely directed vectors in Fig. ^ and are also of the 
order B. In the generalized type II class, B is negative; mi and m2 correspond to vectors of 
order \B\ which are oppositely directed to each other, while rris ranges from very small to 
order \B\. 

Given the bound Eq. (5.9) from the non-observation of {(3(3)oui there are important 
limitations on the allowed values for type I solutions. In this case B is positive so the 
equation 

[m^f = {m^f -B = {m^f - (1.11 eV)^ (5.10) 

will only allow solutions for rris > 1.11 eV. This range is barely compatible with Eq. (5.9). 
Thus the type II case where B < and mi ^ 1712 > seems most probable. 

As an example of a solution for complex M^, consider choosing = 0.2 eV. Then a 
solution is obtained with (compare with Fig. |l|) 

mi ^ ma ^ 1.128 eV,m3 = 0.2 eV,pi ^ 95.1°, /^s ^ 264.9°. (5.11) 

The matrix U, which diagonalizes M,y is obtained from Eq. (5.1), with the approximation 
ujo {(j) = 1, and with now: 

cj-i(r) = diag(0.976 + 0.216i, 0.301 - 0.953z, 0.500 + 0.866i) (5.12) 

This factor introduces CP violation in the {(3(3)o„ process but not in ordinary neutrino 
oscillations. 

Finally the underlying neutrino mass matrix, UMU'^ is "reconstructed" as (in units of 
eV): 

= Mj = 

11 



0.049 - 0.085i -0.855 -0.48n 0.459 + 0.287i 
0.076 + 0.009i -0.025 + 0.131i 
-0.125 + 0.077Z 



(5.13) 



This is structurally similar to the real type II solution displayed in Eq. (4.14), although the 
suppression of the (11) element is not so pronounced. Notice that is considerably smaller 
than the almost degenerate pair mi and Furthermore mi and m2 are large enough to 
possibly play some role in astrophysics. 



VI. SUMMARY AND DISCUSSION 



We investigated the Ansatz Tr {Mi,) = for the underlying (pre-diagonal) three gener- 
ation neutrino mass matrix. It was motivated by noting that in an SO(IO) grand unified 
model where M^, was taken to be real (CP conserving), it corresponds to the well known 
unification of b quark and r lepton masses. While not very predictive by itself it yields 
information complementary to what would be gotten from a complete three flavor analysis 
of all lepton number conserving neutrino oscillation experiments. Specifically from the spec- 
ification of the magnitudes of two neutrino squared mass differences and also of the leptonic 
mixing matrix we can, with some assumptions, find the neutrino masses themselves and 
"reconstruct" M^,. This determination can be sharpened by consideration of the constraints 
imposed by non-observation of neutrinoless double beta decay. 

For the purpose of testing our Ansatz we employed the results of a reasonable best fit to 
all present neutrino experiments (including LSND) by Ohlsson and Snellmann |3^. This fit 
will inevitably be improved in the next few years as new experiments are completed. They 
were able to fit the data without assuming any CP violation. This agrees with assuming 
Ml, to be real. We found two essentially different solutions in that case. The first features 
two neutrinos having approximately equal mass 0.608 eV and a third neutrino of mass 1.217 
eV. This solution is on the borderline of being ruled out by non-observation of (/^/Sqi/)- The 
second solution has two neutrinos with approximately degenerate mass 1.054 eV and a third 
neutrino with a mass 1.36 x 10~^ eV. This solution is very safe from being ruled out by {(3(3qi) 
experiments. It also features a reconstructed My which has an extremely small (11) element. 
Note that, for both solutions, even though My is real there are some (CP conserving) factors 
of i in the mixing matrix when all masses are taken to be positive. Alternatively one may 
have no i's in the mixing matrix while allowing some masses to be negative. The latter form 
is useful for seeing intuitively how Tr(M,^)=0 is possible. 

The case of matching the above best fit data to a complex My was also considered. In this 
situation there are CP violating phases in the lepton mixing matrix which affect the {(3(3)oy 
process but do not affect ordinary total lepton number conserving neutrino oscillations. 
Such phases could also be measurable in principle with the observation of a decay like r~ — »• 
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TT~TT~^~^. The case of complex allows a larger number of solutions. With a simplifying 
assumption there is a one parameter family of allowed neutrino mass sets. Roughly, these 
fall into one of the two types already encountered for real M^. 

A question of some interest is whether the neutrinos are massive enough to play a role 
in cosmology. The relevant criterion I^S] for this to occur is usually stated as X^a^^a > 6 eV. 



For the type II solutions with complex M^, we have found the largest mass sum to be about 
4.5 eV corresponding to mi ^ m2 = 1.62 eV and m-^ ^ 1.22 eV. However this is on the very 
border of acceptability for non observation of {j3l3)ou- 

Future best fits to the neutrino oscillation data can easily be accommodated in the 
present framework. Of course, the predictions for neutrino masses and mixings will depend 
on this input. 
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